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The fate of the Hindu Kush Himalayan glaciers has been a topic of heated debate due to 
their rapid melting and retreat. The underlying reason for the debate is the lack of 
systematic large-scale observations of the extent of glaciers in the region owing to the 
high altitude, remoteness of the terrain, and extreme climatic conditions. Here we 
present a remote sensing—based comprehensive assessment of the current status and 
observed changes in the glacier extent of the Hindu Kush Himalayas. It reveals highly 
heterogeneous, yet undeniable impacts of climate change. 
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Introduction 


The Hindu Kush Himalayan (HKH) region is the freshwater tower of South Asia. It has the 
highest concentration of snow and glaciers outside the polar region, and thus has been 
called the Third Pole (Dyhrenfurth, 1955). The snow and glacier meltwater plays a pivotal 
role in the water supply for those river basins that are arid downstream, in particular where 
there are large irrigation systems that depend on upstream water resources. Mountain 
ranges are particularly sensitive to climate change, and the HKH region is no exception. 
Changes in the glaciers may have a significant impact on the quantity and timing of water 
availability (Rabatel et al., 2013). A comprehensive understanding of the extent and nature 
of changes in glaciers will support downstream hydrological planning and water resource 
management (Rankl, Kienholz, & Braun, 2014). We focus on the development of the 
glacier area over decades, to understand future change and offer a possible evaluation of 
future water quantity and availability. 

There is a large variation in glacier response in the region owing to the diverse terrain 
and climates. There are transitions from colder to drier and a wetter (monsoon dominated) 
to warmer climate from west to east and from north to south. Overall, the glaciers in the 
region are retreating or shrinking and thinning due to climate change (Fujita & Nuimura, 
2011; Kargel, Cogley, Leonard, Haritashya, & Byers, 2011; Kaser, Cogley, Dyurgerov, 
Meier, & Ohmura, 2006; Scherler, Bookhagen, & Strecker, 2011; Zemp, Hoelzle, & 
Haeberli, 2009). Some studies report a different behaviour of glaciers in the Karakoram 
(Hewitt, 2005, 2011; Scherler et al., 2011), which may be explained by a positive trend in 
winter precipitation and limited warming during the melt season (Archer & Fowler, 2004; 
Tahir, Chevallier, Arnaud, & Ahmad, 2011). Climate change also leads to the expansion 
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and formation of glacial lakes, which could lead to an increase in glacial lake outburst 
floods. A number of such floods have already been reported in the region (Bajracharya & 
Mool, 2010; Bajracharya, Mool, & Shrestha, 2007; Richardson & Reynolds, 2000). 

Despite the importance of the HKH region, there is a lack of consistent and 
homogeneous glacier data. Most studies focus on small areas and present a regional 
overview. Several interlinked global glacier inventory initiatives exist, such as the World 
Glacier Monitoring Service (Haeberli, Bosch, Sherler, @strem, & Wallén, 1989), Global 
Land and Ice Measurements from Space (GLIMS) (Raup et al., 2007), the GlobGlacier 
project (Paul et al., 2010), and the Randolph Glacier Inventory (Pfeffer et al., 2014). 
However, none of these initiatives has resulted in a consistent and complete glacier 
inventory for the HKH region. Such a regional glacier database, based on a standardized 
method, is therefore urgently needed as a benchmark against which changes can be 
assessed. Recently, an inventory of glaciers, including debris-covered (DC) glaciers, has 
been developed for the entire HKH region using remote-sensing data (Bajracharya & 
Shrestha, 2011; Bolch et al., 2012). The glacier outlines were compiled using a semi- 
automated remote-sensing method with topographic information obtained from the void- 
filled version of the Shuttle Radar Topography Mission digital elevation model (SRTM 
DEM). The comprehensive assessment of glaciers is based on the status as of 2005 + 3 
years, whereas the extent of glacier change is based on a temporal series of Landsat images 
since the 1980s. 


Study area 


The HKH region extends from 15.95° to 39.31° N latitude and 60.85° to 105.04° E 
longitude, encompassing an area of about 4.2 million km? and extending across all or part 
of the eight countries of Afghanistan, Bangladesh, Bhutan, China, India, Myanmar, Nepal 
and Pakistan (Bajracharya & Shrestha, 2011). There are high concentrations of snow and 
glaciers in the mountains of the region. Meltwater from snow and glaciers feeds the 10 
largest river systems in Asia: the Amu Darya, Indus, Ganges, Brahmaputra, Irrawaddy, 
Salween, Mekong, Yangtze, Yellow and Tarim (Immerzeel, Van Beek, & Bierkens, 2010; 
Kaser, GroBhauser, & Marzeion, 2010). 

Four representative basins were selected for decadal glacier change analysis: the 
Wakhan Corridor in Afghanistan, the Shyok Basin in Pakistan, the Imja Valley in Nepal, 
and the Lunana region in Bhutan (Figure 2). 


Data and methods 
Data 


This study used more than 200 orthorectified Landsat 7 ETM + L1G scenes from the 
Global Land Cover Facility (http://www.landcover.org) from 2005 + 3 years. Images 
with a minimum of cloud and snow cover were selected. The Landsat 7 ETM + images 
suffer from a scan-line corrector failure, causing a wedged-shaped scan-to-scan gap. Prior 
to analysis, the scan-line gaps were filled using a specifically developed software 
extension in the ENVI image-processing package. We used a 90 m resolution SRTM DEM 
to derive the attribute data for each glacier. 

For the decadal glacier change analysis, Landsat 5 MSS, Landsat 7, and Landsat 7 
ETM + images were used. To delineate the glacier outline, the images should have low 
snow and be cloud-free; however, the freely downloadable images rarely meet the ideal 
requirements, hence images must be selected from a range of years rather than a single 
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year. Landsat 5 MSS images from 1976 to 1979 were used to approximate 1980 (hereafter 
referred to as ~ 1980); Landsat 7 TM images from 1988 to 1992, for 1990; Landsat 7 
ETM + from 1999 to 2001, for 2000; and Landsat 7 ETM + from 2009 to 2011, for 2010. 
In each case, the images with the least snow cover and no cloud cover were selected. 


Glacier outline for 2005 


Most glaciers are clearly visible as a clean-ice (CI) surface, but some have debris cover 
(DC) at their tongues. The spectral uniqueness of CI glaciers in the visible and near- 
infrared bands of the electromagnetic spectrum allows the use of simple algorithms. 
However, the delineation of DC glaciers poses challenges because of the non-unique 
spectral signatures and illumination errors due to effects from surrounding materials. 
The images were first segmented using object-based image classification in eCognition 
software (Figure 1) with different parameter settings (type of variable and threshold value) 
for CI and DC glaciers (Bajracharya, Maharjan, & Shrestha, 2014). The Normalized 


LANDSAT and SRTM 


Object-based image classification 









Slope angle 









Classify image objects 










NDVI 8 NDVI 
7) 
3 gg 8 
S Land & Water Mask 9 3 Land & Water Mask ° 
o Ee 5 
ry ® o oO 
2 Slope angle £9 NDSI 5 
DE 
£- 
g 





Area > 0.02 sq.km 


Attribute generation 


Glacier Map/Database 


Figure 1. Flow chart of methodology for delineation of clean-ice (CI) and debris-covered (DC) 
glaciers. Note. NDSI = Normalized Difference Snow Index. NDVI= Normalized Difference 
Vegetation Index. 
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Difference Snow Index (NDSI) (Hall, Riggs, & Salomonson, 1995; Racoviteanu, 
Williams, & Barry, 2008) was used to make an initial selection of segments for CI glaciers. 
In the next filtering step, different variables such as the Normalized Difference Vegetation 
Index (NDVI) (for vegetation), Land and Water Mask (for open water), mean slope and 
mean altitude were determined to identify misclassified image objects, which were 
excluded. The DC glaciers were captured from the segmented image objects using a slope 
threshold value of less than 25°, which captures the entire DC glacier area including 
vegetation, water, snow and surroundings. As DC glaciers only occur below 6000 mas}, 
image objects outside the elevation range of 3000-6000 masl were excluded. Further 
threshold values of NDVI (> 0.3), LWM (50—115.8) and NDSI (=0.005) were used to 
exclude vegetation, snow, and land and water bodies, respectively. The filtering steps were 
not necessarily done in this order. 

These processes satisfactorily delineated the outlines of CI and DC glaciers, which were 
then merged toa single file. Image objects with a surface area of less than 0.02 km? were then 
omitted and the edges of the polygons smoothed. The glacier image objects were transferred 
to a GIS-compatible format and finalized by manual editing at a scale of 1:20,000 by draping 
over respective images and verifying with high-resolution images from Google Earth in a 
GIS environment. Special attention was given to defining the outline and snout of the 
glaciers. CI glaciers require minor manual corrections; DC glaciers need careful manual 
corrections because they closely resemble moraines and bare rock outcrops. 

For the time-series analysis, glacier polygons were generated for 2010 for the Wakhan, 
Shyok, Imja and Lunana areas using the method described above. The glacier outlines for 
the other decades were then generated by manual editing on the recent glacier outline, 
overlaying the Landsat images from 2000, 1990, and ~ 1980. 


Glacier attributes 


Attribute data are assigned to each glacier polygon following international conventions set 
out by the World Glacier Monitoring Service, the GlobGlacier initiative, and the GLIMS 
community (Miiller, Caflisch, & Miiller, 1977; Paul et al., 2010; Raup et al., 2007). For 
each glacier, a GLIMS ID, location, area, elevation range, aspect, average slope, mean 
glacier thickness, and estimated ice reserves were stored in the database. Further details of 
the method are given by Bajracharya and Shrestha (2011). 


Uncertainty 


The uncertainty of the glacier area depends on image resolution and snow cover. 
To minimize the uncertainty, images with the least snow cover were used, and the 
automatically derived glacier polygons were refined manually by draping over high- 
resolution Google Earth images. The maximum offset of the boundary cannot be greater 
than half of the image resolution (i.e. + 15m for TM and ETM + and +30m for MSS). 
The steps to define the uncertainty of the glacier polygon were as follows: 


Buffering of mapped glacier polygons with half of the image pixels 

Calculation of the total number of pixels bounded by each buffered glacier polygon 
Calculation of buffered area by multiplying total pixel count and the area of one 
pixel 

Uncertainty of the mapped glacier taken as the variance in glacier area — that is, the 
difference between actual glacier area (mapped glacier polygon) and buffered 
glacier polygon area 
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The uncertainty of the glacier areas in the present study was 7% for ~ 1980 and 6% for 
1990, 2000 and 2010. 


Results 
Glaciers in the HKH region 


In total, 54,252 glaciers were identified within the HKH region, with a total area of 
60,054 km? and estimated ice reserves of 6127 km? (Figure 2, Table 1). This reveals that 
only 1.4% of the HKH region is glaciated; the total ice reserves are roughly equal to three 
times the annual precipitation (Bookhagen & Burbank, 2006; Immerzeel et al., 2010). 
There is a large variation between river basins; the largest total glacier areas are found in 
the Indus, Brahmaputra and Ganges basins, respectively (Figure 3). 

Over 60% of the total glacier area of the HKH is located in the elevation range 5000-— 
6000 masl (Figure 4). The glaciers below 5700 masl are particularly sensitive to climate 
change unless they are covered by thick debris (Bajracharya et al., 2014). The Indus, 
Ganges and Brahmaputra basins have 79%, 60% and 77% of their total glacier area, 
respectively, below this critical elevation. CI glaciers at low altitude and small glaciers are 
the most sensitive glaciers to climate change in the HKH region. 

A thick debris layer has a strong insulating effect; sub-debris melt rates can be a factor 
of 5 to 10 lower than for CI glaciers (Hagg, Mayer, Lambrecht, & Helm, 2008; Mihalcea 
et al., 2006). A total of 32,000 km” of the glacier area was categorized as either DC or CI 
(glaciers within China were not differentiated). Of this, 9.7% overall was classified as DC, 
and 9.3%, 11.1% and 12.6% of the differentiated area in the Indus, Brahmaputra and 
Ganges basins, respectively. DC glaciers are mostly found at the frontal part of valley 
glaciers and have an average slope of 12° — much less steep than CI glaciers, which have 
an average slope of 25° (Bajracharya & Shrestha, 2011; Bolch et al., 2012). 


HEE Glaciers 
{___] Basin boundary 
[7] HKH boundary 





Figure 2. The Hindu Kush Himalayas (brown polygons), major river basins (black polygons) and 
glacier outlines (blue shaded polygons); red letters indicate the representative glaciers for which 
decadal changes are shown. Note. A = Wakhan Corridor, Afghanistan. B = Shyok Basin, Pakistan. 
C = Inja Valley, Nepal. D = Lunana area, Bhutan. 
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Figure 3. Glacier hypsographs of the major basins in the HKH region. 
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Figure 4. Glacier area distribution in 500 m elevation bands in the HKH region. 


Glacier change in the selected basins 


The time series analysis showed that the glacier area in the Wakhan Corridor is relatively 
stable (Figure 5, Table 2). Some glaciers showed small losses of area of less than 5% in the 
period 1980-1990; there were more noticeable area losses of 4—36% in 1990-2000, 
while most areas were stable or even increased in the period 2000-2010. Over the 30 
years, most glaciers showed a slight reduction in glacier area, which was most prominent 
at the glacier tongues (Figure 5, panel A). Among the individual glaciers, no. 8 increased 
slightly in area over the 30 years; nos. 1 and 6 showed small losses of around 5% of area; 
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Figure 5. Representative examples of glacier change in different parts of the Hindu Kush 
Himalayan region. Note. A = Wakhan Corridor, Afghanistan. B = Shyok Basin, Pakistan. C = Imja 
Valley, Nepal. D = Lunana area, Bhutan. 


the smallest glacier, no. 2, had a marked loss of 27% of area; and all other glaciers showed 
losses of 10-12% of area. 

Most of the glaciers in the Shyok Basin lost area and retreated in the period 1980-— 
1990, showed little change in 1990—2000, and remained static or advanced in 2000—2010. 
Over the 30 years, most glaciers showed some loss of area, but there was almost no change 
in the area of the largest glacier, no. 5. 

The glaciers in the Imja Valley in Nepal and the Lunana region in Bhutan showed 
similar trends, with the highest losses of glacier area observed in the periods 1980-1990 
and 2000-2010. The glaciers in the Himalaya are retreating faster than those in the Hindu 
Kush and Karakorum, with losses over 30 years of 5—55% of area. The West Lhotse 
Glacier in the Imja Valley was the fastest-retreating glacier of all those studied. 


Discussion 


The evaluation of glacier status using remote sensing provides baseline information on 
extension, aspect, slope, elevation, and type of glaciers. Complete inventories for different 
time series are needed to fully understand the likely future challenges for glacier meltwater 
resources and glacial hazards. The glaciers of the HKH region are situated above 3200 
masl, and field-based mapping is difficult given the high altitude, remoteness, harsh 
climatic conditions and lack of logistical support. This limits our understanding of the 
spatial patterns of glacier dynamics and their sensitivity to climate variability across the 
Himalayas at a large scale. Until the launch of Landsat TM in 1982, compilation of glacier 
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inventories was challenging and mostly based on existing topographic maps. Over the past 
several years, remote sensing—based inventories for several parts of the HKH region have 
been entered into the GLIMS database (Raup et al., 2007). However, different methods 
and data sources have been used (topographic maps based on aerial photos with limited 
field inputs, satellite images), with acquisition dates that are years to decades apart. This 
underlines the need for a homogeneous method from a single source of data using recent 
imagery acquired within a short time span. 

The changes in glaciers can be assessed from time series of images. The selected 
basins cover a representative range of climates and glacier systems (Figure 2). A time 
series analysis was performed using a homogeneous source of images (Figure 5, Table 2) 
to assess the glacier change. 

A slight reduction in glacier extent was observed consistently in all glacier tongues in 
the Wakhan Corridor (Figure 5, panel A). The glaciers in the Shyok Basin in the 
Karakorum showed a different pattern, with glaciers advancing and tributaries merging, 
especially during the last decade, consistent with earlier reports (Archer & Fowler, 2004; 
Hewitt, 2005, 2011; Scherler et al., 2011; Tahir et al., 2011). Glaciers in the Himalayas 
(Figure 5, panels C and D) showed consistent retreat in both Nepal and Bhutan. Some ice- 
apron glaciers have disappeared from the mountain slopes, and glaciers are increasingly 
fragmented. These glaciers are all characterized by synchronous ablation and 
accumulation during the monsoon season and are particularly sensitive to climate 
perturbations (Immerzeel, Van Beek, Konz, Shrestha, & Bierkens, 2012). The Imja Lake 
in Nepal (Figure 5, panel C) is one of the fastest-growing lakes in the Himalayas, and the 
DC glacier that connects to the lake has retreated much faster than other nearby DC 
glaciers (Bajracharya et al., 2007). The same phenomenon was observed in the DC 
glaciers that connect to the Raphstreng, Thorthormi and Luggye glacial lakes in the 
Lunana region in Bhutan (Figure 5, panel D). 

The distribution of individual glaciers is important, as small glaciers (Glacier 2 in the 
Wakhan Corridor, Glacier 10 in the Shyok Basin, West Lhotse Glacier in the Imja Valley, 
and Rapthstreng Glacier in the Lunana area) are more sensitive to climate change than 
large glaciers (Bahr, Meier, & Peckham, 1997). A large proportion of the total ice area in 
the Ganges and Brahmaputra Basins is contributed by small glaciers. The distribution of 
the total ice area over different elevation bands is another crucial determinant for climate 
change sensitivity. 

We suggest that in the Hindu Kush glaciers are retreating slowly but steadily, while in 
the Karakoram glaciers show a mixed response, varying from thinning of large DC 
tongues to advancing and surging tributaries. In the Himalayas, glaciers are retreating 
rapidly, with a marked catalytic role played by glacial lakes. However, a complete regional 
understanding will require multi-temporal comprehensive glacier inventories in 
combination with high-resolution climate datasets. 


Conclusions 


A large step has been taken in compiling a reference cryospheric dataset for the entire 
HKH region which is publicly available and can be used as a benchmark against which 
changes can be assessed. The challenge now is to maintain this dataset and to further 
extend it by increasing the temporal interval as far back into the past as available imagery 
allows. 

Larger glaciers (exceeding 15 km* in Shyok and the Wakhan Corridor have shown no 
significant change in area; however, the smaller glaciers are retreating faster. The glaciers 
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in the Wakhan Corridor retreated faster in 1990—2000; the glaciers in the other regions 
retreated faster in 1980—1990. The retreat rate was slower in 2000—2010 in all basins, and 
some glaciers in the Wakhan Corridor and Shyok Basin showed some advance. 

The cryosphere in the HKH region is changing; however, the downstream impact 
remains uncertain. We formulate four key challenges that need to be addressed before 
we can arrive at detailed spatially explicit projections of water availability. (1) 
Precipitation drives the accumulation of glaciers, but spatio-temporal quantifications of 
the present climate, let alone future projections, are highly uncertain. (2) There are no 
fully distributed hydrological models available at the regional scale that include the 
cryosphere with a sufficient level of detail to assess future change in water availability. 
(3) Glacial melt and snowmelt are complex energy-balance processes which are 
influenced by debris cover, black carbon deposition, albedo feedback, and sublimation; 
these processes are poorly understood, and the data needed to derive a full energy 
balance and resulting melt rates is generally lacking at high altitude. (4) There is a need 
for robust algorithms that quantify future glacier extent as a function of model forcing 
and surface features. 
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